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Abstract 

By assuming that not only counter-ions but DNA molecules as well 
are thermally distributed according to a Boltzmann law, we propose 
a modified Poisson-Boltzmann equation at the classical level as start- 
ing point to compute the effects of quantum fluctuations of the elec- 
tric field on the interaction among DNA-cation complexes. The latter 
are modeled here as infinite one-dimensional wires (5-functions) . Our 
goal is to single out such quantum- vacuum-driven interaction from the 
counterion-induced and water-related interactions. We obtain a uni- 
versal, frustration-free Casimir-like (codimcnsion 2) interaction that 
extensive numerical analysis show to be a good candidate to explain 
the formation and stability of DNA aggregates. Such Casimir energy 
is computed for a variety of configurations of up to 19 DNA strands in 
a hexagonal array. It is found to be strongly many-body. 
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1 Introduction 



Dilute solutions of DNA upon the addition of specific multivalent cations 
have remarkable properties pQ. Initially the DNA molecules in aqueous 
solution ionize with two negatively charged phosphates strands remaining, 
giving rise to a highly charged anion with lineal charge density of one neg- 
ative charge per 1.7 A. The system, in the presence of cations, binds them 
by a process known as the Oosawa-Manning (OM) condensation [2] [3] (for 
a review see, e.g., [H HI [5]). When about 90 per cent of the DNA negative 
charge is screened and when the cations have a specific valency +k, usu- 
all}U k = 3 or k = 4, the DNA strands collapse to form rod-like, spheroidal 
and toroidal aggregates, whose size can be experimentally controlled within 
limits [6]. 

There have been different approaches to study these features of DNA 
collapse and structure formation and much progress has been made over the 
last decade. Nonetheless, many aspects remain unclear. Broadly speaking, 
we can identify two avenues that have been explored: 

In one the details of the helicoidal charge distribution of the DNA are 
used to compute the electrostatic interaction between two DNA strands us- 
ing a linearized Poisson-Boltzmann (PB) (or Debye-Hiickel (DH)) approach 
[7] (see also the review article [8j). Within this approach the attraction be- 
tween two like-sign charged but suitably oriented DNA molecules and the 
specificity of cations driving the attraction can be predicted. However, it was 
also realized that when the molecules are three or more frustration comes 
about (for an assembly of several DNA molecules in an hexagonal array on 
a triangular lattice treated within this model see, e.g., [TO] and also |llj). 
Thus, clearly, this force alone cannot account for the formation and stability 
of DNA aggregates. 

In alternative approaches the surface of a single DNA molecule is treated 
as a two-dimensional complex system and statistical mechanical arguments 
lead to the counterion-mediated attraction between two DNA molecules [12} 
[P3"j [T4"l [To] . The key idea there is that the attraction is triggered by local 
correlations and thermal fluctuations not accounted for in the mean-field PB 
approach. These fluctuations being those of the number or charge density 
of the ions, they are constrained to be classical in nature. 

In a recent paper [16] we explored the possibility for quantum vacuum 
fluctuations to be a viable candidate to drive the formation of DNA aggre- 
gates and to hold ("glue") them together. In this paper we want to make 
explicit that analysis and shall present more details and numerical results 
to strengthen that hypothesis. 

1 Counter-ions such as Spermidine 4+ , are surely seen to drive the attraction for DNA. 
Less clear is the matter for k = 2, where, e.g., ions Mg 2+ although they OM bind to the 
DNA strands, are in theory not supposed to trigger attraction [311], while in experiments 
they appear to do that [5]. We shall not address this specificity here. 
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Our starting point is the modification of the PB classical equation to 
include the DNA macromolecules in the Boltzmann thermal distribution. 
This way delta-function potentials appear naturally in the equation to take 
into account the charge distribution of the DNA. Since classical electrostatic 
calculations are available and detailed [7] and since we do not expect our 
one-dimensional (delta- function) model to improve those calculations, we 
focus on the quantum corrections considered here via an effective action ap- 
proach. When the electrostatic potential is small, as the OM condensation 
has taken place, we approximate our modified PB equation to obtain a mod- 
ified DH equation where, as said, the DNA molecule-counterions complexes 
(which, for brevity, we shall often simply call "DNA strands") are modeled 
as delta function potentials (carrying information relative to the charge dis- 
tribution of the strand) as to be expected by general features of self-adjoint 
extension of differential (Hamiltonian) operators [18]. The resulting classical 
action is then taken as the start of a Renormalization Group (RG) type of 
analysis in which time-dependent quantum fluctuations of the electric field, 
propagating with the speed of light in the medium (and representing short 
distance effects) are averaged to generate an effective potential. 

The nice outcome of this procedure is that it produces precisely the 
Casimir energy (for a review see, e.g., [TS], [20] and [5T]) in codimension 2 
(lines in three dimensions) that was obtained in [22J, as we demonstrate. 
There are two scales in the model: the original DH mass (inverse length) 
scale fi, fixed by the parameters of the system, and a new mass (inverse 
length) scale M, introduced through a process of coupling constant renor- 
malization, which is indeed a free parameter. These scales control the range 
of the interaction. 

This energy is calculated for an assembly of N strands (N = 2, 4, 7, 19) 
in arrays in a triangular lattice and we found that: i) it is attractive, irre- 
spective of the charge of the DNA strands; ii) it has a range of O(10) A in 
the simplest model considered; iii) the energy scale in the range is greater 
than ksT; iv) there are important many-body effects. We then conclude 
that it is plausible that this quantum relativistic force is the "glue" holding 
together the aggregate of DNA strands. 

The model presented here has limits. For instance, it suggests that 
the DNA strands collapse to a configuration of zero separation, due to an 
infinitely strong (singular) attractive energy at very short distances, an in- 
stance not occurring in real aggregates. This is an indication that this force 
is only one part of the puzzle and the full picture needs to include mor^E 
It seems to us that the key ingredients to have the full picture are three 
mechanisms: (i) the zero-point quantum interaction, that gives the univer- 
sal attraction ("glue"); (ii) the "frustrated force", that takes into account 

2 One thing to consider within the model is that at very short distances nonlinear 
modifications to the Casimir energy have to play a role. 
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the detailed structure (finite size and helical architecture) and the active 
role of counterions (these alone seem to reproduce the helical architecture 
[17J); (hi) water-related forces, that need to be included. 

In Section [2] we recall the basics of the derivation of the PB equation and 
then propose our modification of it to move on in Section 3 to the quantum 
calculation based on the averaging-over-fluctuations method. The latter is 
done in some details as this should make manifest what sort of interaction 
we are considering here and how we derive it. In Section 0] we present our 
numerical analysis of the interaction between several strands, starting with 
two, and discuss the outcomes and limitations of the model. The last Section 
is dedicated to our conclusions. 



2 The modified Poisson-Boltzmann equation 

Consider the electrostatic potential <J?(cc) due to a charge density at finite 
temperature p(x, T) placed in a medium with dielectric constant e. This 
obeys the Poisson equation 

4-7T 

V 2 $(x) = -— p(x,T) . (2.1) 

Consider now the DNA molecule as a negatively charged rod immersed in 
water at room temperature (T ~ 300K) with dissolved salt whose ions have 
valency z = ±k, with k = 1,2,... (for instance, when the dissolved salt is 
NaCl, k = 1, corresponding to Na + and Cl _ ). The charge distribution of 
the composite system DNA-salt is 

p(x,T) = pbna(x,T) + ke(n + (x,T) - n-(x,T)) , (2.2) 

where pbna(x) is the charge density of the macroion, and we take 

n ± (x,T) = n exv(r^^) , (2.3) 

i.e. the concentration (density) of ions follows a Boltzmann distribution. 
Inserting (|2.2|) and (|2.3p into (|2.ip gives the PB equation 

V 2 cj>(x) = —ken, sinh ( ( 2 . 4 ) 
e V k B T J 

where only the region outside the surface of the macroion is considered, hence 
Pdna = (soon we shall re-introduce pdna)- Equation (|2.4|) is the widely 
used outcome of the mean-field theory that has been extensively applied to 
the study DNA molecules immersed into aqueous media with different salts 
and salt concentrations. Note that both coions (negative) and counterions 
(positive) are present due to the dissolved salt and that in experiments the 
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counterions usually do not come only from the dissolved salt but are added 
separately. 

Among the successes of the PB equation is the prediction of a phase tran- 
sition - governed by the Manning parameter £ = ls/b, with Is = q 2 /eksT 
the Bjerrum length (at room temperature Ib — 7.1 A) and 1/6 the lineal 
charge density that for DNA is b = 1.7 A - where the counterions condense 
onto the DNA strands screening its negative charge to a large extent: when 
£ > 1 counterions stick to the charged rod to form a DNA-counterions com- 
plex. This is the OM condensation [21 E]- 

Our concern is to study the interaction among DNA strands after at 
least 90 per cent of the negative charge has been screened via the OM 
condensation, as this is the reported critical value for collapse. It is then 
reasonable to consider $ small. Furthermore, we explicitly consider the 
charge distribution /)dna and demand that it obeys a Boltzmann distribution 
law as for the ions (see Eq. (|2.3|) ) 



Pdna{x, T) = -71d N a(^)M ex P ( ^T-^r ) > ( 2 - 5 ) 



k B T 

where q < is the charge of the DNA strand with 

N 

u dna(x) = ^2 n{zi)5 2 {x L - k) . (2.6) 

i=l 

This charge density function defines our approximations: we model the DNA 
strands as infinite lines all parallel to the z-axis and located at li in the 
x — y plane with the coefficients Vi(zi) carrying information on the charge 
structure of the DNA strand. We further simplify our model by taking 
v A z i) = v = constant, Vi = 1, ...,N. 

We now put everything together, expand the exponential and stop at 
first order to obtain 



N 

-d 2 z - Vi + /i 2 + a s{2) @i - h 

i=l 



$(£) = J , (2.7) 
which is a modified DH equation. Here fj? = k 2 K 2 , with 



the Debye screening length, 



ek B T 



3 That is why, sometimes, the PB equation above does not have the symmetric expres- 
sion on the right side, sinh(fc), but rather an unbalanced form exp(fci) — exp( — kz)- We 
shall work with the symmetric form (12. 4p . 
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and 

J = -^E^-/-). (2-10) 

6 8=1 

In the following the operator in square brackets will often be indicated for 
brevity as [—V 2 + V(x±)], with obvious notations. 

Delta function potentials are not new in physics (see, e.g., [23]). Their 
appearance in the present case can also be seen as the effect of removing 
points (the locations of the DNA strands) from the domain of the differential 
operator —V 2 . This naturally leads to delta functions as a compensation 
for the self-adjoint extension of such operator [15] . 

Eq. (|2.7|) can be used to compute classical electrostatic interactions. Nonethe- 
less, we shall not do so as such calculations are available and detailed [7] 
and we do not expect our one-dimensional model to improve them. Thus, 
in the following Section, we shall focus on the quantum corrections. 



3 Quantum fluctuations of the electric field 

Let us consider small time-dependent fluctuations: $(x) —* $>(x) + (ft(x,t) 
satisfies the modified DH equation (|2.7p that descends from the action 



a($) = yVzQ$[-v 2 + y(z±)]$ + j$j , (3.ii) 

where we use units fi = c = 1, with c the velocity of light in the medium and, 
for the sake of clarity, we included an integration over time Jq dt even though 
the functions are time-independent. We then demand that to the fluctuation 
field (ft as well is associated an action that is a suitable modification of (|3,lip . 
namely 

A(<ft) = J d A x^cft(-d? - V 2 + V(x ± ))<f> . (3.12) 

Note that in A((ft) the term with the coupling to the "external current" J is 
zero because t 

Jd A xJ<p = J d 3 xJ J dt(f> = 0, (3.13) 

as we impose Jq dtcft = as required for fluctuating fields. The field (ft, 
though, is a quantum field, hence the field configurations that satisfy the 
classical equations (the ones descending from 5A(<ft) = 0) are just on the 
same footing as all other field configurations [22]. The way to consider the 
effects of (ft is to average these fluctuations out to obtain an effective action 
A e f{(&). This is done by considering the generating functional 

Z[$,cft] = f[D<f>]e iA W J[D(ft]e iAw (3.14) 

[D<i>]e~^^' )+corrections ^ , (3.15) 
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where we Wick rotate on the time direction t — > it, and identify A e ff(§) = 
A($>) + corrections. 

Thus we need to compute I = /[D^Je - - 4 ^ that, using standard Gaus- 
sian functional integrals methods (see, e.g., [25]), gives I = [det(— d 2 — V 2 + 
V(x ± ))]-V 2 or 

/ = e X p(-iTrln(-5 i 2 -^ 2 -Vi + y(x ± ))) , (3.16) 

where the identity exp((— 1/2) lndet A) = exp((— l/2)Tr In A) was used and, 
as customary, the determinant and trace have to be computed in terms of 
the eigenvalues of the operator^. To find such eigenvalues we suppose that 

1 f+oo 

cP(x ± ,z,t) = — / dpe^M*^ , (3-17) 
and that (— Vj_ + V(x±))cp p (x±) = E(p p (x±). With this we obtain 

1 = exp {-\ £T ti £T t C dEp{E) Huj2 +p2+ e) ) • (3 - i8) 

Let us focus on I p (E) = (l/2vr) /+~ du \n(u 2 + E 2 ), where E 2 = p 2 + E. 
One has that dI p (E)/dE = l/{2E p ) which gives 

E at (t?i\ i rE i 



llGIlCG 

I = exp (-\ £j ^ jf + °° dEp(E)^E + p 2 ^ . (3.20) 

Thus the corrections to the classical action A(&) are £t with the energy £ 
given by 

1 r+oo A n r+oo , 

£ = 2 J-oo to Jo dE P( E W E + P 2 ■ ( 3 - 21 ) 

This energy is of the form £ = (1/2) "^i i- e - the zero point Casimir energy 
of the system with J2 replaced by (l/27r) J dp j dEp(E) and to by yE + p^. 
Noticing that the density of states p(E) contains information on the location 
li of the N strands present in the system (see later discussion) , clearly d£ jdl 
is a (Casimir) force. 

The mathematical problem of determining the £ of Eq. (|3.2ip has been 
solved in the context of scalar quantum fluctuations for interacting strings 
|22j . Let us give here a brief account of that derivation (see also [26J ) . 



Here a squared length L , which we set to 1, is understood to make the argument of 
the "In" and of "det" dimensionless. Only at the end of the computation (see Eg. (|3. 400 1 
we shall take that into account and shall introduce the proper scale. 
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We first need to find the density of states p{E). Defin^l Hq = — V 2 + p 2 
and H = H + \J2$ {2) {x ~ k) = H + W. The Green's functions are such 
that 

(H-E)G(x,y) =6&>{x-y) , (3.22) 

where < x|G|y >= G(x,y) and similarly for Go(x,y) with Hq. Hence, 
formally we have G = 1/(H — E) and Go = 1/(Hq — E) which allows to 
set-up the Lipmann-Schwinger equation 

1111 

G = — - — = W — - — (3.23) 

H — E Hq — E Hq — E H — E 

= Gq- GqWG = Gq- GqWGq + GqWGqWG = ■■■ (3.24) 

= g »- g » w TVgqW G ^ (3 - 25) 

or 

m \ n i \ , Go(x,li)Go(lj,y) . . 

G{x,y) = G (x,y) + ^ — J , (3.26) 

where a = —A -1 . Now define the energy eigenfunctions as ip n { x ) =< x\E n > 
that give 

oo „ 

p(E) = W d 2 x\Mx)\ 2 8 {2 \E - E n ) , (3.27) 

n=0 

and use 

1 00 1 

G(x, y ) = <*\ H _ E+ie \y>= E M*) En _ E + fe <(y) 

= [V——-+m6(E-E n ))r n (y)- (3-28) 

Clearly 



p(E) = -Im / d 2 xG{x,x) . (3.29) 

7T J 

We need now to prove an identity 

d 2 x < x\— -\k >< -\x > 

Hq — tj Hq — hi 

=< y w^w v ' > - ^ < y^i i < > =^ G »( , "«-( 3 - 30 » 

With this and by using the expression (13.261) for the propagator in the ex- 
pression (|3.29p for the density of states we obtain 



5 Not to clutter the formulae, for the time being we shall not use any special symbol 
for vectors as done earlier. 
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Figure 1: The configuration used for 19 DNA strands where x is the distance 
between nearest neighbors. 

1 8 N 

= Po(£)--^£lnri;- (3-31) 

where po(E) = (l/7r)Im / d 2 xGo(x, x) and Tij = a5ij — Go(k,lj)- By drop- 
ping the /-independent po(E) and by using the fact that is diagonalizable 
we obtain 

1 9 m , = Id 

7T OE 7T oE 

The solution for Go when the energy is Wick rotated to the real 

negative values is [22J 



P( E ) = -Z7^ Trln r 4J = --— lndetf- . (3.32) 



1 



Go(k, lj) = i^ K H E + > (3-33) 

where i^o(a;) is the modified Bessel function of the second kind of order zero 
and lij = |Zj— Hence Gq{L) — > oo when L — > 0. One cures this divergence 
by splitting Gq(L) into a finite part and a divergent part and renormalizing 
the coupling constant A. This process introduces, for dimensional reasons, 
the scale M that, for stability needs to be constrained: M < p |22j. The 
result of this procedure is 



*-G a (L) = 4 + ±l n( ML) + il n (^±Z) (3.34) 

(3.35) 
(3.36) 



where we used the fact that Kq(x) ~ — In a; for small x. We now require 
A — > — oo when L — > so to have a finite A/j = A(L)|i_>o and redefine the 
scale M — > Me 2l / Afl = M which we still require to satisfy M < p. 
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A=1 




5 10 15 20 25 

X(A°) 



Figure 2: Energy of interaction of two DNA strands. The lower (upper) 
curve corresponds to a = 2 (a = 1). Distances x are measured in units 
of 1/fj, ~ O(10) A, while lineal energy density units are estimated to be 
5 x 1(T 2 eV/A. 



Thus the diagonal terms of Tj 7 - are 



fll = h ln \M)^ f Vi = l,-..,N, (3.37) 

and represent the self-interaction. They can be obtained also by taking the 
asymptotic form of Tij 

?ij(kj -> oo) = T5ij . (3.38) 

These terms do not contribute to the Casimir force as they do not depend 
on lij, hence we drop them by considering = Tij/T. Collecting all this 
we can write 

1 /-+ 00 f+oo , Q 

£ = -7^2 y_ oo d p j dE \/ E + p 2 dE ln ( det r «) ' ( 3 - 39 ) 

that when we integrate out the p (using dimensional regularization) and 
partial integrate over E eventually gives 

1 r°° 

£ = — dEIn(detr y ) , (3.40) 

07T Jo 

with 



3 ]n(VE + (J, 2 /M) K J ' 
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Figure 3: Determinant for the two-strand case with a = 2. 



4 Numerical Study of the Interaction Energy 

Our strategy is to study the energy of configurations of DNA strands that 
capture as much as possible the symmetry of arrangements encountered in 
real aggregates [6]. The first case considered is of course that of the two- 
strand interaction. When then proceed to the many-body interactions with 
the strands sitting at the sites of hexagonal lattices like that of Fig. [T] which 
is for the maximum number of strands we were able to consider, i.e. 19 
strands. We also present here results for four DNA strands sitting at the 
vertices of a rhombus. Thus the x-dependance of £(x) is that of the energy 
on the lattice spacing and making x bigger or smaller means to expand 
or shrink the aggregate size, respectively. Of course, for big enough x the 
situation one is describing is that of a dilute solution of DNA molecules, 
which is our starting point. The idealization here consist in the demanding 
a symmetry of arrangement even in the dilute phase. 

To render the expression (I3.40p suitable for such a study we first nu- 
merically perform the integral over E and then plot the resulting expression 
£(x) where 

£( x ) = _Li n ( det7 ( x )) , (4.42) 

and 

*w - S "-Wm il - S,i> <443) 

= Sij - aK (cijx)(l - 6ij) . (4.44) 

The relative distances ly = CijX are expressed in terms of the basic lattice 
distance x and the numerical coefficients qj take the symmetry of the given 
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1.5 




Figure 4: Determinant for the 19-strand case with a = 2. Note that only 
the value x ~ 22 A is the one to consider as the singular point. 



arrangement into account. In (|4.44p we make explicit the choice \i = 1 (i.e. 
the distances are measured in units of The other scale M is constrained 

to be positive and less than fx{= 1) and we write it as < M = e -1 /" < 1. 
In this fashion the range of £ scales with a and we present here results for 
a = 1 and a = 2. 

The two-strand interaction energy is shown in Fig. [2J It is clearly at- 
tractive and finite-range. Similar attractive behaviors for the two-body in- 
teraction have been found in various models [13] . What we observe here 
is that in those models it is not clear why the interaction still needs to be 
attractive for more than two strands. For the Casimir energy we shall soon 
see that this is indeed the case, since this attraction mechanism does not 
suffer of any frustration. 

To establish whether the magnitude and range of this attractive energy 
is indeed relevant for the case of DNA aggregates we need to move, as said, 
to the many-body case. Before doing so we need to first consider that the 
distances are measured in units of ~ O(10) A. This gives a range of 
attraction that can be adjusted by fixing the free parameter a to fit the 
typical distances reached within the aggregates. For the hexagonally packed 
toroidal condensates such distances range between [6j 18 A and 28 A, values 
clearly compatible with the range we obtain for the many-body interaction 
(see later). 

To compare this quantum energy with thermal energy we write the in- 
tegral in Eq. (l3,40p in a dimensionless fashion so that the factor in front is 
ftc/i 2 ~ 5 x 10 eV/A, where c ~ 10 8 ms~ 1 . Hence, as the values of the 
integrals we obtain are 0(1O -3 ), the unit for the energy per length (in A) 
is estimated to be £ ~ 5 x 10~ 2 eV/A. Thermal fluctuations, as computed 



11 




c 
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10 15 20 25 30 35 

X(A°) 



Figure 5: Interaction energy for the four-strand-rhombic case with a = 1 
(upper curve) and a = 2 (lower curve). Distances x are measured in units 
of 1/n ~ O(10) A, while lineal energy density units are estimated to be 
5 x 1(T 2 eV/A. 



for instance in [T3], give for the two-body interaction a maximum value 
of 5 ~ 5 x 10~ 3 eV/A at a separation distance of 10 A. At this distance 
our interaction for the two-body case gives 5 x 10~ 2 eV/A (for a = 1) and 
2 x 10" 1 eV/A (for a = 2), i.e. a result that is between one and two orders 
of magnitude stronger. For the many-body case, the case of importance 
for the aggregates, this factor grows enormously but we cannot trust our 
approximations for x too close to the singular value of the logarithm, say x. 
It is clear, though, that at the distances of relevance this quantum energy 
is stronger (or much stronger) than thermal energy. 

Let us look more closely to the singularity of the energy. We are able to 
numerically evaluate this x for the various cases by plotting the determinant 
that of course shows no such singularity as can be seen in Fig. (x ~ 8 A) 
and Fig. 0] (x ~ 22A) for the two-strand case and for the 19-strand case of 
Fig. [H respectively. That singularity means that if only the Casimir force 
were present the strands would collapse to zero separation, an instance that 
does not occur in the real case because of other forces (not considered here) 
such as the electrostatic force that for more than two strands will give a 
net repulsive effect. Another important factor at such short distances is of 
course the finite size of DNA strands that have a transverse length (radius 
of the cylinder) of lOA. We take the distance x as the limit of validity of our 
approximations. 

Having established the above we have computed the energy for several 
interacting strands having various configurations. We present in Fig. and 
in FigEJthe results for four strands sitting at the vertices of a rhombus and 
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40 45 



Figure 6: Interaction energy for the hexagonal lattice 19-strand configura- 
tion of Fig. [T] with a = 1 (upper curve) and a = 2 (lower curve). Distances 
x are measured in units of ~ O(10) A, while lineal energy density units 
are estimated to be 5 x 10~ 2 eV/A. 



for 19 strands arranged as in Fig. [U respectively. Comparing these plots 
with that of the two strands interaction we clearly see that the attraction 
becomes stronger and acts on a larger range when the number of strands 
increases. 

In the real case of aggregates it is always several DNA strands that 
interact, the two strands being only an idealization. Thus the fact that for 
19 strands we find that (for a = 2) the range of the force is in agreement 
with the typical values reported for DNA aggregates [6j we take it as an 
indication of the validity of our hypothesis that the quantum Casimir energy 
could hold together the aggregates. Furthermore, this force is many-body 
in nature and the many-body effects are big, another reason for taking the 
two-body interaction only as an indication of the real phenomenon. 

That the many body effects are strong we proved in our numerical cal- 
culations where we compared the iV-body energy of Eq. (|4.42|) with that 
obtained by summing up {N/2)(N — 1) two-body interactions. The results 
for four and seven strands are shown in Fig. [7] and Fig. [8j respectively, and 
they indicate that the effect grows with N. 

5 Conclusions 

The main result of this paper is to hint at quantum relativistic effects as 
viable candidates for the collapse of DNA strands into aggregates (after the 
OM condensation has taken place) and for holding them together into sta- 
ble condensates. Rather than approaching the problem by setting-up from 
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Many Body 

Sum of the Parts 




12 14 16 X(A°) 20 22 24 



Figure 7: The lower curve is the interaction energy for the four-strand- 
rhombic configuration (many-body). The upper curve is what is obtained 
by summing-up the 6 two-body interactions. In both cases a = 2. Distances 
x are measured in units of ~ O(10) A, while lineal energy density units 
are estimated to be 5 x 1CP 2 eV/A. 

scratch a model based on quantum electrodynamics we built-up an effective 
model based on the classical PB equation of electrostatics and included the 
DNA strands (modeled here as infinite lines or delta functions) in the Boltz- 
mann distribution, an instance that lead to a modified PB equation. The 
time-dependent fluctuations we then studied are quantum in nature, propa- 
gate at the speed of light in the medium and give rise to a Casimir force that 
we studied in various settings. In particular, we focused our attention on the 
difficult problem of computing such interaction for several DNA strands and 
were able to overcome the analytical challenges with numerical calculations 
performed for a variety of cases, most of which with hexagonal symmetry of 
arrangement (the typical situation reported in experiments). 

The numerical calculations show many interesting features: the interac- 
tion is attractive and short range for the two-body case and is many-body, 
the departures from the "sum of two-bodies" being important and grow- 
ing with the number of strands; the distance at which our approximations 
stop working is also obtained as the value at which the determinant func- 
tion vanishes; finally, the magnitude and range of the interaction is such 
that it could explain the formation and stability of DNA aggregates, as a 
preliminary comparison with reported data shows. 

Our model is a primitive one and to make full contact with experiments 
we propose it as one part of the puzzle as it needs to be seen as one of the 
concurrence of three mechanisms: (i) the zero-point quantum interaction, 
that gives the universal attraction ( "glue" ) ; (ii) the "frustrated force" , that 
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Figure 8: The lower curve is the interaction energy for the 7-strand configu- 
ration, i.e. for 7 strands sitting at the vertices and at the certer of a regular 
hexagon (many-body). The upper curve is what is obtained by summing-up 
the 21 two-body interactions. In both cases a = 2. Distances x are measured 
in units of l//i ~ 0(10) A, while lineal energy density units are estimated 
to be 5 x 1(T 2 eV/A. 

takes into account the detailed structure (finite size and helical architecture) 
and the active role of counterions (these alone seem to reproduce the helical 
architecture [17]; (iii) water-related forces, that need to be included. That 
the interaction (I3,40p alone does not give the full picture is seen from the 
singularity at a value x such that the determinant becomes zero. If only this 
force were present the DNA-cation complexes would collapse to zero sepa- 
ration, but at very short distances we should include nonlinear corrections 
to this force itself, and, even leaving aside water-related forces, at short dis- 
tances the Wigner crystal effects should become important. Our goal here 
was to single-out the role of the zero-point quantum vacuum energy and we 
have shown that it is quite plausible that it plays a key role in the onset of 
the formation and in the follow-up stability of DNA aggregates. 

It is pleasant to see that quantum effects might be essential for under- 
standing an important biological problem (other quantum effects are im- 
portant for enzyme catalysis [27] or speculated to be important for neural 
activity [28]) as this might serve as a solid basis for a more general under- 
standing of the role of quantum mechanics for life [29] . 
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